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Abstract 

Quark number susceptibilities at finite quark chemical potential are inves- 
tigated in the framework of the Polyakov-loop-extended Nambu Jona-Lasinio 
(PNJL) model. A detailed comparison is performed between the available 
lattice data, extrapolated using a Taylor expansion around vanishing chem- 
ical potential, and PNJL results consistently obtained from a Taylor series 
truncated at the same order. The validity of the Taylor expansion is then 
examined through a comparison between the full and truncated PNJL model 
calculations. 

Investigations of QCD thermodynamics using lattice simulations predict a tran- 
sition from the hadronic to the quark-gluon phase around a critical temperature 
of 0.2 GeV, at vanishing quark chemical potential \x q . This transition is signalled 
by a steep rise in energy density, pressure and entropy density as functions of the 
temperature, and known to be a crossover [1] when quarks are included. 

Based on model calculations [2] as well as lattice QCD simulations [3], the exis- 
tence of a critical endpoint in the phase diagram is suggested. The precise location 
of this point is still under debate. It can possibily be identified by observables that 
are sensitive to singular parts of the free energy [4], such as quark number suscep- 
tibilities. Such susceptibilities have been studied both in lattice calculations [5,6], 
and using phenomenological models [7-10]. Both diagonal and off-diagonal suscepti- 
bilities can serve as useful diagnostics to elucidate the nature of strongly interacting 
matter [11]. 
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The first lattice results for quark number susceptibilities at finite chemical poten- 
tial have been obtained in [5] through the Taylor expansion method. An expansion 
of the pressure in powers of \i q jT around /i q — was performed up to fourth order 
in n q /T. Quark number susceptibilities were then obtained as second derivatives 
of the pressure with respect to the chemical potential. The presence of a peak in 
these observables at large chemical potentials has been interpreted gnal of 

growing fluctuations in the baryon density, as the critical endpoint is approached. 
Improvements were implemented in [6] by expanding the pressure up to sixth order 
in n q /T. In this case, the peak was shifted to a smaller temperature, and a dip 
was observed at T ~ 1.05 T c which, together with the increased error bars in the 
extrapolated lattice data, made the presence of the peak less convincing. 

In the present paper we investigate quark number susceptibilities within the 
Polyakov-loop-extended Nambu Jona-Lasinio (PNJL) model [12-15] which has re- 
cently been used successfully in comparisons with a variety of lattice QCD data 
[16-21]. In this model, quarks propagate in a temporal background gauge field rep- 
resenting Polyakov loop dynamics, while developing at the same time a dynamical 
mass through their coupling to the chiral condensate. Here we present results for 
the quark number susceptibilities at finite chemical potential. A comparison is first 
performed with truncated Taylor expansions in \i q jT derived from lattice QCD. 
Then the convergence properties of the Taylor series are examined by comparison 
with the full, non-truncated PNJL result. 

The Euclidean action of the two-flavor PNJL model is 



where ip is the Nf = 2 doublet quark field, a = 70 7 and 74 = 270 in terms of the stan- 
dard Dirac 7 matrices and mo = diag(m u ,m,d) is the quark mass matrix. V(ip,ip ') 
is an 577(2) x SU(2) invariant four-fermion interaction acting in the pseudoscalar- 
iso vector /scalar-isoscalar quark-antiquark channel 1 



The quarks move in a background color gauge field = A 4 = iA , where A = 
S^o gA^ t a with the SU (3) c gauge fields A^ and the generators t a = \ a /2. The matrix 

lr Thc scalar diquark interaction that has been considered in [20] is not relevant in the present 
discussion since the range of chemical potentials explored here is lower than the critical /i q for 
diquark condensation. 




Here 7i is the fermionic Hamiltonian density given by: 



U = -i^ (a ■ V + 74 m - 0) t\) + VO, ^) , 



(2) 




(3) 



2 



valued, constant field relates to the (traced) Polyakov loop as follows: 



, /3 

Pexp 



(ij drA^ ^Tre^, (4) 



In a convenient gauge (the so-called Polyakov gauge), one can choose a diagonal 
representation for the matrix 0, 

= 03 A 3 + 08 A 8 , (5) 

which leaves only two independent variables, 3 and 8 . 

The term —(V/T) U of the action (1) involves the effective potential W($, $*, T) 
which controls the thermodynamics of the Polyakov loop. An improved form for this 
effective potential was introduced in our previous work [20]. It involves the logarithm 
of «/($), the Jacobi determinant which results from integrating out six non-diagonal 
SU(3) generators while keeping the two diagonal ones, 3j8 , to represent $: 

^( $ '^*' T ) = --a(T) $*$ + b(T) In [1 - 6 <F$ + 4 ($* 3 + $ 3 ) - 3 (<T$) 2 ] (6) 
J. 2 

with 

a(T) =a + ai (^j + a 2 (^j * , b(T) = b 3 (^j ' . (7) 

The logarithmic divergence of W(<£>, $*,T) as $, $* — > 1 automatically limits the 
Polyakov loop $ to be always smaller than 1, reaching this value asymptotically 
only as T — > oo. The parameters a« and 63 are taken from Ref. [20], where they were 
determined to reproduce lattice data for the thermodynamics of pure gauge lattice 
QCD up to about twice the critical temperature 2 . The resulting parameters are 

a = 3.51 , ax = -2.47 , a 2 = 15.22 , b 3 = -1.75 . 

The critical temperature T for deconfinement in the pure gauge sector is fixed at 
270 MeV in agreement with lattice results. 

The NJL part of the model involves three parameters: the bare quark mass which 
we take equal for u- and c?-quarks, the coupling strength G and a three-momentum 
cutoff A. We take those from Ref. [14] : 

m u4 = 5.5 MeV , G — 10.1 GeV" 2 , A = 0.65 GeV , 

fixed to reproduce the pion mass and decay constant in vacuum and the chiral 
condensate as = 139.3 MeV, / w = 92.3 MeV and (4> u ^ u ) = -(251 MeV) 3 . 



2 At much higher temperatures, where transverse gluons begin to dominate, the PNJL model is 
not supposed to be applicable. 
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After performing a bosonization of the PNJL action and introducing scalar and 
pseudoscalar auxiliary fields, a and 7?, the PNJL thermodynamic potential becomes: 



2 

n (T, N , a, $, <F) = U ($, $*, T) + ^ (8) 

_2AT f / _^_{TlnM +3 $e- ( ^-^ )/T +3$*e- 2( ^-^ )/T +e- 3( ^-^ )/T l 
J 7 (2tt) 3 l l j 

+ T In [1 + 3 $* e -^ + ^ T + 3 $ e - 2 ^ + ^ T + e" 3( ^ + ^ )/T ] + 3 AE P 9 (A 2 - p 2 ) } 

where the quark quasiparticle energy is E p = \/p 2 + m 2 and the dynamical (con- 
stituent) quark mass is the same as in the standard NJL model: m = mo — a = 
itlq — G{^ip). The last term in the previous equation involves the difference AE P 
between the quasiparticle energy E p and the energy of free fermions (quarks). It is 
understood that for three- momenta \p\ above the cutoff A where NJL interactions 
are "turned off' , a is set to zero. 

In general, the Euclidean action Se is formally complex in the presence of the 
temporal gauge field (j) [22,23]. It is real only at vanishing chemical potential, 
fi q = 0. At nonvanishing chemical potential one has in addition (<&) ^ (<&*) [24]. 
Fluctuations beyond mean field are at the origin of ($) ^ (<&*) for /j, q ^ [25] (see 
the more detailed discussion in [20]). This paper deals with self-consistent solutions 
and predictions of the mean-field equations 

with if = a, 03, 8 . At this mean field level, the additional constraint of 0j being real 
implies that the action is minimized by 8 = 0. It follows that $ ~ Tr exp(iA 3 <fiz/T) 
is real (i.e. $ = $* at the mean field level). The temperature dependence of the 
Polyakov loop and chiral condensate, obtained from Eq. (9) can be found in Ref. [20]. 

A Taylor expansion of the pressure in fi q /T has been performed in Ref. [20] in 
order to compare with corresponding lattice data [6]. The pressure is written as a 
series 



rp4 



n=0 



CO = 



1 d n (p(T,fi q )/T 4 



9 {nJTy 



(10) 

n q =o 



The Taylor expansion coefficients C2, C4 and cq, calculated in the PNJL model, turn 
out to agree very well with the corresponding lattice data. 

The PNJL model calculations can provide useful insights concerning the con- 
vergence of the expansion (10). The sign problem at non-zero chemical potential 
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restricts lattice QCD thermodynamics to extrapolations including only the first few 
terms of Eq. (10), whereas the model is not restricted to such truncations. In 
Ref. [18] we have shown that, for the pressure, very good agreement is reached be- 
tween the full result and the truncated one, already at fourth order in n q /T and up 
to fi q /T ~ 1. In the case of the quark number density, the agreement is still very 
good at small chemical potentials, but discrepancies are observed at large fJ> q /T in 
the vicinity of the phase transition. 

In the present paper we investigate the quark number susceptibilities defined as 

The isospin-symmetric case is considered, with equal chemical potential for up and 
down quarks. The full result for Xq is compared with the truncated one given by 
the Taylor expansion 

= 2c 2 + 12 c 4 (^) 2 + 30 c 6 (|)\ (12) 

In Fig. 1 we show a comparison between truncated PNJL result obtained from 
Eq. (12), and the corresponding lattice data, consistently expanded to the same 
order in fi q /T. One finds remarkably good agreement between the lattice data and 
the PNJL predictions. At large chemical potential [i q /T ~ 1 we observe, consistently 
with the lattice data, a pronounced peak sligthly below T c , followed by a dip around 
T ~ 1.05 T c . This feature was also discussed in [6]. 

Fig. 2 presents a study of the convergence of the Taylor expansion in ft = \i q jT . 
A comparison is made between full PNJL results and those truncated at fourth and 
sixth order in ft, for a sequence of chemical potentials 0.4 < fi q /T < 1. For low 
chemical potentials, the agreement between the full result and the truncated one is 
rather good. Only in the case of [i q jT = 0.6 and in the region of T around the phase 
transition, the first discrepancies start to appear. 

The right panel of Fig. 2 shows a comparison between the full PNJL result (con- 
tinuous line), the one truncated at fourth order (dashed line) and the one truncated 
at sixth order (dotted line) for fi q /T = 1. Evidently, the Taylor series converges 
badly in the region around the phase transition. The full result shows a peak, not 
a divergence, around T/T c ~ 0.9, indicating that the critical temperature for the 
phase transition decreases with increasing chemical potential. The coefficient c 6 is 
sizeable around the phase transition, so that the sixth order contribution differs 
significantly from that of fourth order. The dip above T c occurring at sixth order 
is entirely an effect of the truncation of the series. It is absent in the full result. 
It reflects the corresponding dip exhibited by cq. in the fourth order result (dashed 
line) it is in fact absent as well. The height of the peak in the full result is much 
reduced as compared to the corresponding one found at sixth order. No divergence 
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Figure 1: Scaled quark number susceptibility \q as function of T/T c for different values of Hq/T. 
The curves are PNJL model results obtained from the Taylor series expansion in n q /T around 
fiq = (see Eq. (12)). Corresponding lattice data are taken from Ref. [6]. 



is observed in the full result. In fact, our model predicts that the phase transition 
is a crossover, for this value of fi q . This analysis indicates that the truncated Xq can 
be misleading in the sense that it incorrectly suggests the possibility of a first-order 
transition which is not realized in the full susceptibility. 

From the present work we can draw the following conclusions. The available 
lattice results for quark number susceptibilities at finite chemical potential, obtained 
from a Taylor expansion in fi q /T around \i q = 0, can be successfully reproduced in 
the framework of the PNJL model through an analogous Taylor expansion truncated 
at the same order in n q /T. The comparison between the full PNJL result at finite 
chemical potential and the truncated one shows good convergence of the Taylor 
series up to fi q /T < 0.4 — 0.6. For larger values of fi q /T, significant discrepancies 
are observed between full and truncated results in the region around the phase 
transition. The singular behaviour of the susceptibilities observed in the sixth-order 
result is not present in the full calculation. The transition predicted by our model 
for /J,q/T = 1 is still a crossover, reflected by the finite height of the peak in the 
corresponding quark number susceptibilities. 
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Figure 2: Left: comparison between the full and truncated PNJL results for two different (small) 
chemical potentials. Right: comparison between full and truncated results, at different orders, for 
fJ, q /T= 1 (sec Eq. (12)). 
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